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Abstract 

We describe a method for analyzing the phase space structures of Hamiltonian 
Q \ systems. This method is based on a time-frequency decomposition of a trajectory 

^ I using wavelets. The ridges of the time-frequency landscape of a trajectory, also 

pH . called instantaneous frequencies, enable us to analyze the phase space structures. 

In particular, this method detects resonance trappings and transitions and allows a 
characterization of the notion of weak and strong chaos. We illustrate the method 
with the trajectories of the standard map and the hydrogen atom in crossed mag- 
. . netic and elliptically polarized microwave fields. 
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1 Introduction 



^ \ The phase space of a typical Hamiltonian system is a mixture of periodic, 

^ ' quasiperiodic and chaotic trajectories, and the stable and unstable manifolds 

of these trajectories. The analysis of the trajectories of those systems is central 
to many branches in physics, ranging from celestial mechanics to chemistry. 
Several methods have been designed to discriminate between regular (periodic 
or quasiperiodic) and chaotic trajectories. For two degrees of freedom, the most 
commonly used method is a Poincare section which is a plane section of the 
three dimensional energy surfaces. These sections give pictures of phase space 
and allow to analyze qualitatively the stability properties of the system by re- 
vealing, e.g., the invariant tori, the stable islands surrounding elliptic periodic 
orbits and the spatial extent of chaotic trajectories. However, these sections 
do not provide very accurate information and the visual inspection of these 
sections for higher dimensional systems is less straightforward since, e.g., for 
Hamiltonian systems with three degrees of freedom, these sections are four 
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dimensional. 

Much of our understanding of nonlinear systems is based on their frequencies 
and especially resonances between them (perturbation theorems, intramolec- 
ular dynamics [1,2], Laskar's work in celestial mechanics [3,4,5] and particle 
accelerators [6]). One way to obtain accurate information on the stability of the 
system is to associate with a given trajectory an indicator such as a frequency 
(in general, several indicators are associated with a given trajectory) and con- 
struct, e.g., a frequency representation of phase space. In this frequency space, 
regular regions are characterized by a very small variation of the frequency 
as the time interval on which this frequency is computed increases (diffusion 
constant in Laskar's work [7]). The main advantage of this method is that it 
is generalizable to systems with any number of degrees of freedom [3]. 
However, a frequency representation of the phase space does not provide the 
history and the local properties of a trajectory. The computation of asymptotic 
quantities (such as Lyapunov exponents, entropy or fractal dimensions) does 
not provide dynamical information. These quantities only reflect the asymp- 
totic behavior of the system. For many purposes infinite-time quantities are 
not sufficient : rather than knowing whether or not a trajectory is chaotic, it 
is naturally desirable to say when, where, and to what degree a trajectory is 
chaotic. 

In this article, we describe a method that gives a time-frequency lepresenteition 
of trajectories, i.e., that provides dynamical information on the system. The 
idea is to follow a given trajectory and see if the phase space structures where 
the trajectory passes by are revealed by analyzing a single trajectory. For 
quasiperiodic and chaotic systems, time-frequency representation is revealing 
because it shows when and how motion can be trapped in some resonance 
zones associated with quasiperiodicity. It also show the degree of chaoticity of 
some part of phase space visited by the trajectory. 

In a nutshell, we associate a set of time-varying frequencies with the trajec- 
tory by decomposing it on a set of elementary functions, the wavelets, which 
are localized in time and frequency. Wavelet decomposition has been used to 
analyze trajectories of Hamiltonian systems in celestial mechanics [8] and in 
molecular dynamics [9,10]. In Refs. [8,10], the main frequency was extracted 
by computing the frequency curve where the modulus of the wavelet transform 
is maximum. However, stopping at this main frequency can give misleading or 
wrong information about resonance transitions. Instead, we have found that 
the relevant information of the trajectory is contained in the ridges of the 
landscape formed by the magnitude of the coefficients of the time-frequency 
decomposition [11]. A ridge curve (which has a certain time length) corre- 
sponds to time-frequency energy localization and is also called instantaneous 
frequency curve. We show that the instantaneous frequencies point to reso- 
nance trappings and transitions, and give a clear characterization of the de- 
gree of chaoticity. 
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In Sec. 2, we describe the time-frequency method and the extraction of ridges 
by two methods: the windowed Fourier transform and the continuous wavelet 
transform. In Sec. 3, we present briefly the two models which are used to il- 
lustrate the method in Sec. 4: the first one is the standard map which is used 
as a benchmark model for the method since the phase space structures of this 
model are well-known, and the second one is the hydrogen atom driven by a 
elliptically polarized microwave field which is a basic and realistic problem in 
atomic physics with three degrees of freedom. 

We show how time-frequency analysis reveals the phase space structures: For 
two degrees of freedom, we compare this analysis with the Poincare surface of 
section, and we show that the time-frequency analysis can be carried out for 
higher dimensional systems. 

In Appendix A, we give an approximation of the ridge curves for a quasiperi- 
odic trajectory and in Appendix B, we compute the effect of finite time on 
the computation of the ridges. 



2 Time-frequency analysis 

Given a real signal f{t) which can be, e.g., one coordinate of the system (one 
position or one momentum for instance), two methods are used to compute a 
time-frequency (denoted u - ^) representation of /: windowed Fourier trans- 
form and continuous wavelet transform. We assume that we have computed 
the trajectory for all times t G M. The effect of finite time interval is ana- 
lyzed in Appendix B. For more details on time-frequency analysis, we refer to 



2.1 Windowed Fourier transform 



The windowed Fourier transform (also called Gabor transform) oi f{t) is given 



where the window g is chosen to be a Gaussian window g{t) — e~*^/^'^^/(cr^7r)^/*^. 
We notice that in Laskar's work [3,4], the Hanning filter x{t) = 1 — cos(27rt/T) 
over a time interval [0, T], is chosen rather than a Gaussian filter. The Hanning 
filter has the advantage of a finite support whereas the Gaussian filter has a 
minimal time-frequency resolution. We will consider an energy density in the 
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time-frequency plane called spectrogram 

Psf{u,0 = \Sfiu,0\'. 

Wc notice that Psf{u,—^) — Psf{u,$,) for / real. It is thus sufficient to 
consider the positive frequency part of the time-frequency plane. 

2.1.1 The algorithm 

The algorithm to compute Sf{u,^) uses one fast Fourier transform (FFT) of 
/(■ + u)g{-) at each time u since 

/+00 
fit + u)g{t)e-'^'dt. 
-oo 



2.1.2 Time-frequency resolution 

The resolution of the transform is constant in time and in frequency. The time 
spread around a point {u,^) in the time-frequency plane, defined as 

+00 \ 1/2 



where gu,(,it) = g(t — u)e~^^^, is independent of u and For a Gaussian window 
g{t) = e"* /^"^ /(ct^tt)^/^, it is equal to (7/-\/2. The frequency spread around 
{u, ^) defined as 

/ I r+00 \ 1/2 

t^/req(M,0 = {uj-0^\gu,i{uj)\'^dwj 

does not depend on u and ^ since gu,s,{^) = — i)e^'^^^~^\ For a Gaussian 
window, it is equal to l/{a^/2). We notice that the product of the time spread 
with the frequency spread at a given point (u, ^) in the time frequency-plane 
is constant and larger than 1/2 (Heisenberg uncertainty). It is minimum for a 
Gaussian window. 



2.2 Continuous wavelet transform 



The continuous wavelet transform of f{t) gives a time-scale representation of 
the trajectory and is given by 

wf{u, = i= £J f{t)r (^) dt, (2) 
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where the mother wavelet ip is chosen to be a Gabor (modulated Gaussian) 
wavelet, also called Morlet-Grossman wavelet: ip{t) = e"'^e^'^ /^'^ /(a^Tr)^/"^. 
The transform depends on rj which is the center frequency of the wavelet 

1 f+OO 



ZTT J-oo 



where il>u,s{t) — s — u)/s]. The time- frequency representation is ob- 

tained by the relation between the scale s and the frequency ^: 

We will consider the normalized scalogram 

Pwf{u,C-v/s) = -\Wf{u,s: 

s 
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which can be interpreted as the energy density in the time-frequency plane. 
Again, we notice that Pwf{u, — ^) = Pwf{u,^). For the computation of the 
wavelet decomposition, we only consider the positive frequency part of the 
time-frequency plane. 



2.2.1 The wavelet algorithm 

In order to use a fast wavelet algorithm, we consider that / and are periodic 
(and hence Wf{u, s) is periodic in the variable u) with period T, the time span 
of the trajectory. By taking the Fourier transform of Eq. (2), one has: 

Since V^(a;) is known exactly, one has to compute / by FFT and then retrieve 
Wf by an inverse FFT. For this computation, the scales are chosen according 
to a dyadic samphng Sij — 2'+^/"T'/A'" for i = 1, . . . , m and j — 1, . . . ,n where 
m is the number of octaves, n is the number of voices per octave and is 
the number of points in the time series on the interval [0,T]. The relative 
error produced by the algorithm on the computation of the frequency is then 
Acu/cu = « (log2)/n. This error decreases when the number of voices 

per octave is increased. We chose n ^ 300 is our computations (which gives a 
relative error of 2 ■ 10~^). In what follows, we choose m — log2 N — 3 where N 
is the total number of points in the series. 

In order to compute the continuous wavelet transform, we use the routines of 
WaveLab, written in Matlab [16]. 

Remark: In practice we choose rj = 2 and a = 6. We notice that the wavelet 
does not satisfy the admissibility condition ^(0) = which is required for 
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reconstructing the trajectory. However ^(0) = (47r(T^)^/^e"'^^''^/^ is negligible 
since at] is large enough (e.g, if ar] = 12, |'?/'(0)| < 10^^°). If one wants to 
satisfy the admissibility condition, the wavelet has to be chosen according to 

with c = (7rcT2)-i/4(i _2e-3'^''?V4 + e-<x2^2^-i/2^ r^j^^ relation between the scale 
and the frequency is changed. It becomes ^ = (1 + e"'^ ^ )v/^ up to the first 
order. For the parameters we choose, this error in determining the frequency 
is neghgible. 



2.2.2 Time-frequency resolution 

The time-frequency resolution depends on the scale of the wavelet. The time 
spread around a point {u, s) in the time-scale plane, defined as 

/ r+oo \ 1/2 

CTtimeiu, S) = U {t - uf\iju,s{t)\'^dt] , 

is independent of u (the moment of observation of the frequency content) and 
is proportional to the scale s = r//^ for a mother wavelet of the form ip{t) — 
e"'*g'(i). For a Gabor wavelet, aume is equal to ar]/{V^^). From this equation, 
we see that for low frequency the time spread of the wavelet is large, whereas 
for high frequencies, it is small. This allows a more accurate computation 
of the frequency since for low frequencies a larger time interval is needed 
(large enough to contain several oscillations at that frequency). Increasing 
the parameters a and rj increases the time spread which varies as 1/^. The 
frequency spread around ^ — r]/s defined as 

/ \ P + OD „ \ 1/2 

afreqiu,s) = J^^ {uj - ri/sf\tlJu,si^)\^dujj , 

is independent of u and is linearly dependent on the frequency for a mother 
wavelet of the form ip^t) = e"'^g{t). For a Gabor wavelet, cr/^eg is equal to 
^ / {ar]\/2) . We notice that the product of the time spread with the frequency 
spread at a given time- frequency point is independent of u and ^ — r]/s and 
is larger than 1/2 (Hciscnbcrg uncertainty). The Heisenberg uncertainty is 
minimum for the Gabor wavelet. Increasing or decreasing the product of the 
parameters err/ is a compromise between the time resolution and the frequency 
resolution. 
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2.3 Differences between windowed Fourier transform and wavelet transform 

Both methods have time-frequency resolution hmitations for the determina- 
tion of the instantaneous frequencies. The windowed Fourier transform relies 
on an a priori choice of length of the window a. Any event (trapping, transition, 
etc.) happening on short time-scales (less than a) or with small frequencies 
(less than a~^) is missed by this method. On the contrary, the main advantage 
of the wavelet method is that it follows the rapid variations of the instanta- 
neous frequencies since it adapts the length of the window according to the 
frequency [12]. 

For quasiperiodic trajectories, there are practically no differences between both 
methods. However, for chaotic trajectories where variations of frequency with 
respect to time are expected, the wavelet basis is more adapted since it leads 
to a better time-frequency resolution. In what follows, we choose the wavelet 
basis in order to compute the time-frequency content of the trajectories. 

2.4 Ridge plots 

The different main frequencies of the trajectory can be obtained by looking at 
the ridges of the spectrogram or normalized scalogram, also called instanta- 
neous frequencies. These ridges are local maxima with respect to the frequency 
^, of the energy density in the time-frequency plane which is the spectrogram 
or the normalized scalogram. 

For a periodic trajectory f{t) = e**^*, the spectrogram 

P5/(H,O = 2<70Fe-<^'(-«)' 
and the normalized scalogram 

are maximum for ^ = uj, independently of the time u. Therefore, the ridge 
plot will present a single fiat ridge located at ^ = uj. 

For a quasi-periodic trajectory f{t) = J2k ^k^^'^''^ , the spectrogram and the 
normalized scalogram present a sum of localized peaks and interference terms : 
The wavelet transform of / is given by the sum of the wavelet transform 
of ^fce"^*'* since the wavelet transform is linear. Therefore, the normalized 
scalogram is a sum of peaks located at C, = ^^k and of width proportional to 
u!k/{T]<^), with some additional interference terms of the form 

The effect of these terms is visible if (a) at least two amplitudes /l^ and Ai are 
large enough, and (b) if the difference of their frequencies uJi—uJk is smaller than 
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uji/{r]a) (see Appendix A). In order to distinguish two frequency components 
uji and cijfe of amplitude Aj^ and Ai of the same order, the parameters of the 
wavelet must be chosen according to the condition r]a > {cui + u!k)/\'-^i — i^k\- 
The time-frequency landscape has approximately a set of flat ridges (if the 
interferences are negligible) at ^ = a;^. 

Remark: Laskar's frequency map analysis. In order to determine accurately 
the frequencies of a quasiperiodic trajectory and the location of the ridges, we 
use Laskar's frequency map analysis [3,7,4] which is a Fourier based analysis of 
quasiperiodic functions f{t) — J^k OfcC*'^''* known on a flnite interval [0, T]. The 
approximation of the main frequency cui , called u![ is obtained by maximizing 
the modulus of 

where xit) — 1 — cos(27rt/T) is a Manning filter. The amplitude a\ of the com- 
ponent with frequency u[ is obtained by projection on e^'^i*. Laskar proved 
that this method gives the main frequency with an accuracy of order [4] 
(i.e., uj[ — uji = 0(T~^)) instead of an ordinary Fourier transform which gives 
an accuracy of order T~^. The approximation of the next frequency 002 is ob- 
tained by iterating the above procedure using /i(t) = f{t) — a'^e^'^i*. Since the 
basis {e^'^fc*} is not orthogonal for the scalar product (•,•), a Gramm-Schmidt 
orthogonalization is necessary to obtain the next frequencies o;^ for k > 3. 



For a chaotic trajectory, a ridge is a curve or a segment of curve ^lodu) in 
the time-frequency plane which is at each time u a local maximum of the 
normalized scalogram, i.e.. 



-^Pwfiu,0 



-0, 



C=?ioc(") 



< 0. 



?=^!oc(w) 



Each ridge has a weight which is the value of the normalized scalogram on this 
ridge (it varies continuously in time) . We will refer this value as the amplitude 
of the ridge. We call main ridge or main frequency, the ridge ^rn{u) (or a set 
of ridges) where the normalized scalogram is maximum: 



Pwf{u,^miu)) 



max Pwf{u,i). 



In what follows, we will detect ridges by analyzing one coordinate of the sys- 
tem. We will show that the set of ridges of one coordinate contains all the 
frequency content of the trajectory, that is to say, the ridge plot will not be 
dependent on the choice of the coordinate used for the computation of the 
time-frequency transform. The respective amplitudes of the ridges depend on 
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the specific coordinate chosen for the computation of the normahzed scalo- 
gram. For instance, the main ridge depends in general on the chosen coordi- 
nate whereas the ridge plot does not. 

Detection of the ridges. For each time u, we are searching for the local maxima 
of the spectrogram or the normalized scalogram. In general, there is an infi- 
nite number of these ridges. We only detect the ones with sufficient amplitude 
according to a threshold e, i.e., we determine the set of ridges ^ioc{u) such that 



3 Models 



The time-frequency analysis can be performed on trajectories of maps or of 
fiows. In this section, we illustrate the time-frequency analysis on the two 
models: the well-known standard map and the hydrogen atom in crossed mag- 
netic and elliptically polarized microwave field. 

The standard map is a two-dimensional area-preserving map {x, y) i— > (x', y') : 



x' = X — a sin y (4) 
y' = x' + y mod 2tt. (5) 

This example provides a benchmark for the time-frequency analysis since the 
phase space structures are well-known in this example. 

The hydrogen atom driven by a elliptically polarized microwave field (polarized 
in the orbital plane) and a constant magnetic field perpendicular to the orbital 
plane is described by the following Hamiltonian [17]: 

Px+P^ 1 

H{x,y,p^,Py,t) = -^—^ ; ,) , ■> + VF{x,y,u;t) + VB{x,y,p^,Py), (6) 

2 +y 

where Vp is the contribution of the elliptically polarized microwave field and 
Vb is the one of the magnetic field : 



VF{x,y,u!t) — F{x cos cut + ay sin cut), 
B 

VB{x,y,p^,Py) = —{xpy-yp^) + — (x^ + 

where uu is the frequency of the microwave field and a is the ellipticity degree. 
For a = 0, the microwave field is linearly polarized (along the x direction), 
and for a = 1, it is circularly polarized. In the circularly polarized case, 
the problem has two effective degrees of freedom. For a g]0, 1[, the field is 
elliptically polarized and is a three degree of freedom problem. 
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The magnetic field B perpendicular to the orbital plane of the electron limits 
the diffusion of trajectories [18,19]. With it in place, trajectories will stay for 
a longer time close to the nucleus before ionizing, allowing longer trajectories 
to analyze by time-frequency analysis. 

Rescaling of the Hamiltonian: By rescaling positions, momenta and time, 
we can assume that uo = 1. The rescaling of time is obtained by consider- 
ing the Hamiltonian uj^^H{x, y,Px,Py, t/ur) where H{x, y,Px,Py, t) is Hamilto- 
nian (6). Then we rescale positions and momenta : x' — lo'^^^x, p'^ — u;~'^^^Px: 
y' — up-l^y and = uj~'^^^py^ which is a canonical transformation generated by 
S{x\ y',Px,Py) = ~^~'^^^{x'px + y'py)- Then we rescale the momenta by chang- 
ing H into uj^^^H(uj~^^^px,uj~^/^py,x,y). We notice that this transformation 
does not change the equations of motion. The resulting Hamiltonian is the 
same as Hamiltonian (6) with cu — 1 and a rescaled value for the amplitude of 
the microwave field and the magnetic field: The amplitude of the microwave 
field is rescaled as F' — Fuj''^^^ and the magnetic field diS B' — B/u. In what 
follows, we assume that uj — 1. 



3.1 Semi-parabolic coordinates 



The singularity of the potential in Hamiltonian (6) is removed by changing 
the coordinates to the semi-parabolic coordinates and a nonlinear change of 
time. We perform the canonical transformation generated by 

S{u,V,Px,Py) = Px - UVPy. 

The expression of the old coordinates as functions of the new coordinates is 



V? — 



X 



2 ' 
y = uv, 

_ UPu - VPv 
Px 9 I 9 ' 

_ VPu + UPy 

Py- ^2 + ^2 ■ 

The Hamiltonian becomes 



H= ^"/^l, - , ^ , +F (— — —cost + auvsint 



B B^ 
+ -iup^-vp^) + —{u^ + v^)^. (7) 
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Wc consider t as a new variable w with conjugate momentum pyj. The new 
Hamiltonian becomes time-independent: 



^ = ^T7^^^ - 2 I 2 +-^ — n — cosw + auvsmw 

B 

We denote E the energy of the system described by H' (which is conserved 
since H' does not depend on time). Then wc change time by multiplying 
the Hamiltonian by v? + v"^. The new time r satisfies dt/dr = v? -\- (sec 
Ref. [20]). We consider the effective Hamiltonian whose equations of motions 
are the same as the ones of Hamiltonian (7): 



2 2 

^e// = ^4^ + (p.-i?)(«^ + ^^)-2 

^F{y? + v^^ \ — ^ — cos w + auv sin w \ 
B B^ 

+ -(1^2 + ^;2)(^ip^ - + _(ti2 + ^2)3^ 

where i^e// = 0. We shift the momentum p^ by a factor —E and the effective 
Hamiltonian becomes 



^e// = ^^+P.(^^^ + ^^)-2 



2 , 2n 



+F(u + V ) y — - — cos w + auv sin w j 

B B^ 
+ -{u^ + v^){up^ - vpu) + —{u^ + vy, (8) 

We perform the following rotation of coordinates: 



_ w w 

u — u COS V V sin — , 

2 2 ' 

. w w 

V — —u sm \-v cos — , 

2 2 



which is a canonical transformation generated by 

S{u,v,w,Pu,Pv,Pw) = -pu[ucos{w/2) - vsm{w/2)] 

—py[usm{w/2) +vcos{w/2)] — wpu 
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The Hamiltonian (8) becomes 



~ 2 I ~ 2 
fV Pu I Pv ,-2 
Heff = — " 



smw -\- uv cos w 



smw 



(9) 



where H^ff is still equal to zero. In the circularly polarized case {a — 1), 
the above Hamiltonian is independent of w, thus is constant, i.e., p^j + 
{upy — vpu)/'2 = const. Thus the circularly polarized problem has two degrees 
of freedom. Poincare surfaces of section will be drawn: the surfaces are the 
ones with Pp = (in the rotating frame), i.e., upu + vp^ = 0. 



3.2 Choice of initial conditions 



We choose the initial conditions for the circularly polarized problem. Given 
the values of F and B, the energy of the system is chosen according to Ref. [18] 
in the maximum configuration. At this energy, there is a stationary point of 
the flow and interesting dynamics happens in this case (see Ref. [18] for more 
details). We notice that the stationary points belong to the surface Pp — 0. 
The maximum configuration is obtained for v — 0, i.e. for 

p~^^^{l-B)4-Ful, (10) 
where Uq is the positive real solution of 

^{i-b)4-f4-a^o. 

We notice that the equation that determines always one positive real 

root and two complex conjugated roots for F > and B < 1. 

Once we have fixed F. B and the energy of the system p^, we compute Poincare 
sections by the following procedure. These surfaces will be in the (,t, t/)-plane. 
We consider one point (x, y) on this plane. We compute the associated u and 
v: 



12 



u — y + 1/^ + x 



1/2 



V = (sgn y) (yx^~+y^ - 



1/2 



The following conditions upu+vpy — (which is Pp — 0) and Heff{u, v,Pu,Pv) 
determine p„ and p^: 



2 



-D \ 

1 - — j {u^ + v^) +£5{u,v) 
B \ 

1 - — j {u^ + v'^) + e5{u,v] 



where 



5{u, v)= (1 - B){u' + v'f - 8p~^ - 4F{u' - v') + 16{u'^ + v^)'^ 



1/2 



and £ = ±1 depending on whether we specify the trajectory crossing the 
Poincare section Pp = in one direction {Pp > 0) or the other (Pp < 0). 
We use Poincare sections in order to exhibit phase space structures and show 
that the time-frequency analysis of chaotic trajectories is able to reveal these 
structures. 



4 Illustration of the time-frequency method 



In what follows, we consider the trajectories of the standard map (discrete 
time) and the trajectories of the hydrogen atom in crossed fields (continuous 
time). We will show that the phase space structures encountered by a tra- 
jectory are revealed by the time-frequency analysis of a single coordinate of 
this trajectory, i.e., by the computation of the ridges of the time-frequency 
landscape, the instantaneous frequencies. 

We compute the time-frequency content of real time series (e.e., one coordi- 
nate of the system) obtained by iterating the map or integrating the flow. 
We only compute the time-frequency content of the trajectory in the positive 
frequency region. 



4-1 Trajectories of the Standard Map 



4.1.1 Quasiperiodic trajectory 

We consider a trajectory of the standard map for K = 0.9710 starting at 
Xq — 4.176562 and yo = over a time span of 512 iterations (which corresponds 
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°0 500 1000 1500^ 2000 2500 3000 

Fig. 1. Phase portrait (first panel), normalized scalogram (second panel) and ridge 
plot (third panel) of a quasiperiodic trajectory of the standard map obtained for 
K = 0.971 and initial conditions Xq = 4.176562 and yo = 0. 



to T = 3217 for the corresponding flow). This trajectory, plotted in the first 
panel of Fig. 1, corresponds to the invariant torus (golden mean torus) whose 
break-up is studied in Ref. [21] using frequency map analysis. The second panel 
of Fig. 1 depicts the normalized scalogram in the time-scale plane, i.e. the 
amplitude of the wavelet coefficients, obtained from the x coordinate. From 
this time-scale (or time-frequency) landscape, we extract the local maxima 
which define the ridges or the instantaneous frequencies of the trajectory. The 
third panel of Fig. 1 shows the different ridges of the normalized scalogram (the 
ridge detection threshold is e = 0.1). These ridges are approximately constant 
in time as one expects. The upper curve represents the main frequency ujq = 
(3 - V5) /2. The other two main frequency curves in the ridge plot are located 
at cui = 1 — 2u!o and UJ2 = 3a;o — 1. For a system with two degrees of freedom, 
there are at most two independent frequencies: in this case, they are equal to 
cuq and 1 which is the natural frequency for constructing the Poincare map 
from the flow (kicked rotor). Each ridge (a priori an infinite number) is a 
combination of these two frequencies. 



We can also see the interferences (oscillations around the expected frequency) 

between these frequencies. In order to reduce these oscillations, we can increase 
the parameter rja. We notice that the amplitudes or frequencies of these oscil- 
lations do not depend on the length of the time interval in which the trajectory 
is computed (see Appendix A). 
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Fig. 2. Phase portrait (first panel) and ridge plot (second panel) of a weakly chaotic 
trajectory of the standard map obtained for if = 1.2 and initial conditions xq = 2.5 
and yo = tt. 

4.1.2 Weakly chaotic motion 

We consider a trajectory obtained for the standard map with K = 1.2 and 
initial conditions xq = 2.5 and yo = tt. We analyze this weakly chaotic tra- 
jectory over the time interval of total length T = 12868. The upper panel 
of Fig. 2 depicts this trajectory. The lower panel depicts the main ridge (the 
ridge detection threshold is e = 0.5) and shows clearly different trappings 
as time evolves: trappings in the resonance 1:3 for 150 < t < 500, in 2:5 for 
5600 < t < 6800, in 4:7 for 8700 < t < 9300, and in 3:5 for 11200 < t < 12600. 
These trappings result from the trajectory passing nearby islands surrounding 
elliptic periodic orbits (the other ridges contain the secondary frequencies of 
the nearby quasiperiodic motion). The transitions between these resonances 
occur quite smoothly, passing by other resonance trappings. They occur when 
the chaotic trajectory passes nearby hyperbolic periodic orbits where great 
variations of the frequencies are expected. The trajectory is weakly chaotic 
since the main information can be obtained from a single instantaneous fre- 
quency curve. In this regime, looking at the frequency where the spectrogram 
or scalogram is maximum (main instantaneous frequency) is meaningful. 

4.1.3 Resonance transitions 

The time-frequency resolution of wavelets allows us to analyze the mechanism 

of resonance transition. Figure 3 shows the normalized scalogram in the time- 
scale plane and the ridge plot (with ridge detection threshold e = 0.1) obtained 
for the X coordinate of a trajectory of the standard map with K = 1.5 and 
initial conditions xq — 1.35 and yo = on a time interval [1700, 2500]. In this 
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Fig. 3. Ridge plot of the x component of the trajectory of the standard map with 
K = 1.5 and initial conditions xq = 1.35 and yo = on the time interval [1700, 2500]. 
The upper panel is the time-scale plane and the lower panel is the ridge plot. The 
bold line denotes the main frequency curve. 



figure, the bold curve represents the main ridge or the main frequency (i.e. 
where the normalized scalogram is maximum, see Sec. 2.4). 



In this time interval, several transitions occur. For instance around t = 1860, 
the main frequency jumps discontinuously from the ridge located around 
cui ~ 0.33 to the ridge located around 002 ~ 0.136. One could wrongly con- 
clude that a sharp resonance transition occurred in the system. By looking 
at secondary ridges, i.e. the ridges with smaller amplitudes than the main 
frequency, one can see that both ridges existed before and after t = 1860; 
Nothing discontinuous happens in the system. Prom the ridge plot, the mech- 
anism of resonance transition appears to be the following one: A second ridge 
at ^ = CJ2 increases in amplitude from t = 1700 tot = 1860 while the first ridge 
at ^ = uji decreases in amplitude. On some time interval (from t = 1800 to 
t = 2000), several ridges of the same order in amplitude appear. At t = 1860, 
the second ridge becomes the main frequency. This situation is typical from 
other ridge plots of trajectories showing transitions from ridges and hence 
resonance transitions (see Fig. 12). Thus wc have a new characterization of 
resonance transition that applies in more general situations. Resonance tran- 
sitions are a manifestation of the fulfillment of the Chirikov resonance overlap 
criterion [22]. 
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Fig. 4. Ridge plot of a strongly chaotic trajectory of the standard map with K = 5 
and initial conditions xq = 0.95 and yo = 0- 

4-i-4 Strongly chaotic motion 

When K increases, the trajectories are increasingly chaotic (filling a bigger 
part of phase space). More and more short ridges appear, and in those cases, 
looking at the main frequency leads to erroneous results as shown above since 
several ridges have the same amplitude. The timc-frcqucncy information can- 
not be reduced to (or deduced from) a single instantaneous frequency. For in- 
stance, the main frequency curve could follow a short fiat ridge and one could 
wrongly deduce a trapping region. We consider a strongly chaotic trajectory 
of the standard map obtained with K — 5 and initial conditions a;o = 0.95 and 
yo = 0. Figure 4 shows the ridge plot obtained from the normalized scalogram 
(with a ridge detection threshold e = 0.5). It shows that no trappings occur 
in this time interval. The presence of a lot (but a finite number) of ridges 
reflects the broad-band spectrum of a chaotic trajectory. These pictures lead 
to a characterization of weak and strong chaos by looking at the number of 
important ridges in the system: Weak chaos is characterized by one main con- 
nected instantaneous frequency curve whereas strong chaos is characterized 
by multiple short ridges (the number of them increases as K increases) . 

4-2 Hydrogen atom in crossed fields 

We consider in this section the trajectories of the Hamiltonian flow of the hy- 
drogen atom driven by crossed magnetic and elliptically polarized microwave 
field in the rotating frame represented by Hamiltonian (9) in the maximum 
configuration. 

Otherwise specified, we consider in this section the following values for the 
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Fig. 5. Poincare section of the Hamiltonian (9) in the maximum configuration for 
F = 0.015 and B = 0.3. 

parameters: F = 0.015 and B = 0.3. Figure 5 depicts a Poincare section of 
the Hamiltonian in the circularly polarized case. We analyze the trajectories 
by applying the wavelet transform on one real coordinate x{t) or y{t). 

4-2.1 Periodic orbit 

We consider a closed trajectory in phase space obtained for initial conditions 
(x, y) (0.4835, 0.5335) in the Poincare section. A projection of this periodic 
orbit on the plane {x, y) is shown in the upper part of Fig. 6. The ridge plot 
of this orbit is depicted in the lower panel of Fig. 6. 

The ridge plot displays only one fiat ridge located at the value of its frequency 
loa ~ 0.2952. The harmonics are not visible on the ridge plot since these ridges 
are too small compared to the main ridge (ridge detection threshold e = 0.1). 
The bumps at m = and u = 4500 are a consequence of the finite time 
interval on which the periodic trajectory is computed. This effect is discussed 
in Appendix B. 

There are other elliptic periodic orbits in the system: For instance, the one 
obtained from the initial conditions in the Poincare section {x, y) — (—1.014, 0) 
with frequency ujb ~ 0.3033, and for (x,y) = (—0.807,0.558) with frequency 
uc ~ 0.2260. 

4-2.2 Quasiperiodic trajectory 

The ridge plot of a quasiperiodic orbit obtained for initial conditions {x, y) = 
(0.7, 0.5) is depicted in Fig. 7. It displays a set of three main ridges around uji ~ 
0.17, a;2~0.25 and ou^^O.^A (with a ridge detection threshold e = 0.1). The 
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Fig. 6. Ridge plot of a periodic orbit of Hamiltonian (9) for F = 0.015 and B = 0.3 in 
the maximum configuration obtained for initial conditions {x^y) = (0.4835,0.5335) 
in the Poincare section. The upper part of the figure represents a projection of the 
periodic orbit on the x — y plane. 
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Fig. 7. Ridge plot of a quasiperiodic orbit of Hamiltonian (9) for F = 0.015 
and B = 0.3 in the maximum configuration obtained for initial conditions 
(x, y) = (0.7,0.5) in the Poincare section. The upper part of the figure represents 
the Poincare section of the trajectory. 



oscillations are due to the interactions between ridges. This effect is explained 
in Appendix A. The frequency map analysis gives the following values for the 
first frequencies: ui — 0.1735, ijJ2 — 0.2531 and — 2uji — 0.3470. 
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Fig. 8. Ridge plot of a chaotic orbit of Hamiltonian (9) for F = 0.015 and B = 0.3 
in the maximum configuration obtained for initial conditions: {x,y) = (—0.5,0.7). 
The upper part of the figure represents the Poincare section of the trajectory. 



4-2.3 Chaotic orbits: resonance trappings and transitions 



We consider a trajectory in a chaotic region obtained for initial conditions 

(x, y) = (—0.5, 0.7) in the Poincare section. Figure 8 depicts a Poincare section 
of this trajectory and the ridge plot of the coordinate x{t) over the time interval 
[0, 2.85 10"^] (with a ridge detection threshold e = 0.5). 

Figure 9 depicts an enlargement of the ridge plot of Fig. 8 on the time interval 
[1.35 10^,2.65 10^] of two coordinates x and y with a smaller ridge detection 
threshold e = 0.1. Similar ridge plots are obtained with the other coordinates 
Px and py. 

Three different type of regimes are observed on these figures: trapping regions 
where the frequency is constant, weakly chaotic regions where the frequency 
is changing smoothly (for instance, from t = 1.6 10"^ to t = 1.95 10^), and 
strongly chaotic regions characterized by many short ridges (for instance be- 
tween t = 5.5 10^ and t = 8.5 10=^). 

The two main trapping regions in these figures are obtained from t = 10^ 
to t = 1.38 10^ where the trajectory is trapped around the elliptic island 
surrounding the periodic orbit with frequency ujb, and from t = 2.23 10'^ to 
t — 2.35 10*^ where the trajectory is trapped around the elliptic island with 
frequency cuc- 

The part of phase space visited by the trajectory is compact (this is due to the 
magnetic field) so we do not expect great variations of the main frequency. 
In Fig. 10, we depict a resonance transition that occurs in the system for 
F — 0.02 and B = 0.3 (still in the maximum configuration) for initial con- 
ditions {x,y) = (—0.7, —0.5) in the Poincare section. This figure displays the 
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Fig. 9. Ridge plot of the chaotic orbit of Fig. 8 on the time interval 
[1.35 10'',2.65 10"^]. The Tipper part of the ridge plot is obtained by analyzing the 
coordinate x{t) and the lower part is the one for the coordinate y{t). 
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Fig. 10. Ridge plot of the chaotic orbit of Hamiltonian (9) in the maximum config- 
uration with the parameters F = 0.02 and B = 0.3 obtained for initial conditions 
{x, y) = (—0.7, —0.5) in the Poincare section. The upper part of the ridge plot is 
obtained by analyzing the coordinate x{t) and the lower part is obtained from the 
coordinate y{t). 

ridges of the coordinates x and y with a ridge detection threshold e = 0.1. 
There are two trapping regions: from t — 3000 tot — 8000 and from t — 10400 
to t = 14500. The first trapping region is characterized by two main frequen- 
cies uJi = 0.1750 and 002 = 2c<Ji, and the second trapping region has three 
main frequencies u![ = 0.4106, Lj'2 = 0.4630 and uj'^ = 2uj2 — t^i- Both regions 
are separated by a strongly chaotic region characterized by many short ridges 
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Fig. 11. Poincare section of tlic trajectory of Fig. 10. The upper left panel is the 
trajectory from t = to t = 6.1 10^, the upper right panel is the trajectory from 
t = 3 10^ to t = 8 10^, and the lower panel is the trajectory from t = 10.4 10^ to 
t = 14.5 10=^. 

(resonance transition region). Figure 11 displays Poincare section of the trajec- 
tory from t — to t — 61400 (upper left panel), the section of this trajectory 
which is trapped between t = 3000 and t = 8000 (upper right panel), and the 
section of the trajectory trapped from t = 10400 to t = 14500 (lower panel). 
In the trapped regions, the trajectory nearly fills a one-dimensional part of 
the Poincare surface. 

From Figs. 9 and 10, we notice that the ridge plots are the same for the x 
coordinate of the y coordinate. The only difference is that a given ridge has 
different amplitudes according to whether it has been computed using x{t) or 
y{t). One exception is Fig. 10 from t — 3000 to t — 8000: only the ridges of 
sufficient amplitude are detected (with a ridge detection threshold e = 0.1) 
and the ridge around ^ ^ 0.35 which is the main ridge for the x coordinate 
has a low amplitude for the y coordinate. This ridge appears if we include low 
amplitude ridges in the detection. 

The noticeable differences between the ridge plots for x and the ones for y are 
in the strong chaotic region where there are no elliptic islands and thus no 
trappings. In this regime which is the same as the one in Fig. 4 for the standard 
map, the significance of a frequency or a set of frequencies as an indicator is 
not clear. It means that the basis is not adapted to the strong chaotic regime. 
However, computing the ridge curves is useful since it provides a way to lo- 
cate these strongly chaotic regions. It also provides a quantitative distinction 
between weak and strong chaos by looking at the number of short ridges. 
We highlight that computing a single instantaneous frequency from the maxi- 
mum wavelet coefficients can provide erroneous information. Figure 12 depicts 
the ridges in a strongly chaotic zone of phase space (the same trajectory as for 
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Fig. 12. Ridge plot of the same trajectory as for Fig. 8 on the time interval 
[1200,2000]. The bold curve represents the main frequency curve. 

Fig. 8 on the time interval [1200, 2000]). The bold curve is the main frequency 
ridge. By looking at this main frequency, we cannot say if the trajectory is in a 
strong chaotic region or if it sharply jumps between different trapping regions. 
The computation of secondary ridges shows that it is indeed a strongly chaotic 
zone without trappings. 

4-2.4 Time-frequency analysis of higher dimensional systems 

The time-frequency analysis can be carried out for higher dimensional sys- 
tems. We consider the elliptically polarized case with a = 0.9 represented by 
Hamiltonian (9) in the rotating frame. In this case, Pw is no longer a conserved 
quantity so Hamiltonian (9) has three degrees of freedom. Figure 13 depicts 
the ridge plot of the x coordinate with a ridge detection threshold e = 0.5, 
and Figure 14 depicts the ridge plot of the y coordinate with a ridge detection 
threshold e = 0.1. 

We notice that from this figure the trajectory looks more chaotic than in 
the circularly polarized case which is in agreement with Refs. [23,24] (see the 
differences between Fig. 8 and Fig. 13). From Figs. 13 and 14, we detect reso- 
nance trappings and weakly chaotic regions. For instance, from t = 1.45 lO"' to 
t = 1.65 10^ the trajectory is trapped in a region where the main frequencies 
are cui — 0.127 (which is the main frequency for y) and UJ2 — 2a;i (which is 
the main frequency for x). The same trapping is observed on a longer time 
interval from t = 2.65 10^ tot = 3.15 10^. This is confirmed by a Fourier anal- 
ysis of the segments of the trajectory: Figure 15 shows the power spectrum 
of two segments of the trajectory: from t = 1.45 10*^ to t = 1.65 10^ which 
is expected to be close to quasiperiodic, and from t — 3.2 10^ to t — 3.6 10^ 
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Fig. 13. Ridge plot of the x coordinate of the chaotic orbit of Hamiltonian (9) with 

OL = 0.9 with the parameters F = 0.015 and B = 0.3 obtained for initial conditions 
{x,y) = (—0.5,0.7) and solution of Eq. (10); Px and py are solution of H^ff = 
and xpx+ypy = 0. The upper part of the ridge plot is for the time interval [0, 2.4 10^] 
and the lower part is for [2.4 10^,4.9 10^]. The ridge detection threshold is e = 0.5. 
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Fig. 14. Ridge plot of the y coordinate of the same trajectory as Fig. 13 with a ridge 
detection threshold e = 0.1. 

which is expected to be strongly chaotic with a broad-band spectrum. 



By increasing the ridge detection threshold (see Fig. 14), we clearly see the 

trapping regions (with one main flat ridge) and the strongly chaotic regions. 
We notice that there are regions where there is one main flat ridge with addi- 
tional short ridges (for instance from t = 10"^ to t = 1.25 10*^ or from t = 3.6 10^ 
to i = 4 10^). A Fourier analysis of these segments of trajectory shows a main 
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Fig. 15. Power spectrum of two segments of the trajectory of Fig 13 (x coordinate): 
from t = 1.45 10^ to t = 1.65 10^ (upper panel), and from t = 3.2 10^ to t = 3.6 10^ 
(lower panel). 

frequency uj ~ 0.15 and a broad-band component. This is a nearly trapped 
region near some elliptic island. Prom Fig. 14, we notice that the trajectory 
is trapped in this island from t = 4.22 10^ to t — 4.3 10^, and the main 
frequencies are cui — 0.15 and UJ2 — 0.22. 



5 Conclusion 



In summary, we showed how instantaneous frequencies based on a ridge ex- 
traction of a wavelet decomposition of a single trajectory coordinate reveal 
the phase space structures (resonance transitions, trappings, etc.) of chaotic 
systems. Our method also gives a quantitative characterization of weak and 
strong chaos through trappings and through the number of short ridges. We 
have showed that this time-frequency analysis can be carried out on trajec- 
tories of maps and flows. Note that nothing in our analysis depends on the 
dimensionality of the system: Given a time series (which can be e.g., one of the 
coordinates of the system whatever the number of such coordinates may be, or 
indeed data sets) , we can compute the ridges of the time- frequency landscape, 
leading to the instantaneous frequencies. Moreover, the computational time is 
independent of the dimensionality of the system. The time-frequency method 
based on a single trajectory analysis is therefore very well suited for explor- 
ing phase space structures of various systems independently of the number of 
degrees of freedom. 
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A Time-frequency analysis of quasiperiodic trajectories 



Given the following quasiperiodic function 

/(t) = e-i* + ^afce-'=*, (A.l) 

k>l 

where e C, we derive an approximation of the ridges of the time-frequency 
landscape obtained using the windowed Fourier transform and the continuous 
wavelet transform. 



A.l Windowed Fourier transform 



The windowed Fourier transform of trajectory (A.l) is (up to an overall con- 
stant (T-\/27r): 

fc>i 

We want to study the behavior of Sf{u, ^) around ^ ~ We expand Sf{u, C.) 
up to second order in the variable x = ^ — cui : 



where 



^y(O) _ gii^iM _|_ ^ Q^giu)fcUg-o-2(u,j._t^i)2/2 
k>l 



k>l 



The ridge around uji is the local maximum of the absolute value of Sf{u,^) 
for each u. The expression of the spectrogram is 
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2Re + ] + 0{x'). 



+ 

The ridge is approximately located at 

Re 



2Re(,S/(o),5/(2)*) + |5/(i)| 



When the terms a^e^'^^'-'^'^^'^^-'^/^ are small (compared to 1) which is the case 
when Gk is small or uj^ is sufficiently far from uji, the expression of the ridge 
becomes : 

i{u) = ^1 + - uJi)Ak cos[(wfe - uji)u + 0fc]e-^'('^'=-'^i)'/2^ (A.2) 

fc>i 

where at = AkC'^'^^. The instantaneous frequency curve is again a quasiperiodic 
function with frequencies Uk — wi. It oscillates around the main frequency 
ui with amphtude proportional to e~"' (^k-^^i) /2 xhese oscillations can be 
reduced by increasing the length of the window of the transform a. 



A.2 Wavelet transform 



The wavelet transform of trajectory (A.l) is (up to an overall constant cr \/27r): 
Wf(u, s) = e^'^^^e-^'^'^'^i-"/^)'/' + afee^'^'="e-^'"'('"^-^/")'/'. 

k>l 

We want to study the behavior of the normalized scalogram Pwf{u, ^) around 
s ~ r]/cAJi. We expand s~^/'^W f {u, s) up to second order in the variable x — 
SLUi/fj — 1 : 

s-^/^Wf{u, s) = + Wf^^^x + Wf^^^x^ + 0{x^), 

where 



k>l 



2 2 

—ar) 



k>i ^1 



1 2 2 f^k , 

1 — a Tj 1 

Va;i 



The ridge around uii is given by the local maximum of Pwf{u, ^) for each u. 
The expression of the normalized scalogram 
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+ [2Re (wf''^Wf<^>) + ] x' + 0{x^). 



The ridge is approximately located at 



Re (]Vf ^^Wf'>) 



2Re(W^/(o)W^/(2)*) + |W^/W|2' 



When the terms a^e 



a2r,2(u;fe/a;i-l)2/2 



are small (compared to 1) which is the 



case when ak is small or uok is sufficiently far from cui, the expression of the 
ridge becomes : 



^{u) = u;i + E M<^k/coi - l)Ak cos[K - + 0fc]e-^''''('^'=/'^i-i)'/2. (A.3) 



As for the windowed Fourier transform, the instantaneous frequency is a 
quasipcriodic function with frequencies ujk — uji. It oscillates around the main 
frequency ui. The amplitude of the oscillations is proportional to e^"' i'^k/i^i-'i-) /2_ 
They can be reduced by increasing the parameters of the wavelet a and rj. 



B Effect of finite time 



In this section, we investigate the effect of the finite time on the computation 
of the instantaneous frequencies. We consider the signal f{t) — e"^* on [0, T]. 



B. 1 Windowed Fourier transform 

The windowed Fourier transform is computed by considering a periodic signal 
/ with period T constructed from /, i.e. f{t) = f{t) for t G [0,T[ and f(t + 
T) = f{t) for all t. However the signal / may not be continuous if the value 
of the original signal / at and T are not the same. The discontinuities lead 
to an alteration of the instantaneous frequencies near these points. 
We compute the ridge near a point of discontinuity, e.g., u — 0. The windowed 
Fourier transform of / is equal to 



k>l 
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For u close to zero, two terms arc predominant in the sum: n = — 1 and n — 0. 
On the interval [— T, 0], the function / is equal to e*'^(*+^). Thus 

Jo Jo 
We introduce the error function erf(2;) = ^ Jq e~^^dt: 



a 71 



1/4 



+e 



iwT 



^ ^ ] gi(u;-«)Wg-c72(a;-0V2 

[el(I±ii±-(^).e.| 



-u — ia^{u — ^) 



We assume that T is large compared to u and (T^|a; — We introduce the 
complex error function w{z) as w{z) — [1 + erf {iz)] (see Ref. [25]). Since 
erf z — > 1 when z — > oo and |arg z\ < tt/A, the expression of the windowed 
Fourier transform becomes: 



w{x + iy) + e^'^w{—x — iy) 



where x = a{u — ^)/\/2 and y = —u/{a^/2). We notice that when cuT = 
mod 27r, the expression of Sf becomes equal to 

Sf{u,0 ~ (47r(72)V4e-(-«)e-^(-0V2^ 

since w{z) + w{—z) — 2e~^^ . That is the windowed Fourier transform of the 
periodic signal e*"^* on R. When uT ^ mod 2tt, the ridge is no longer 
constant and located at ,^ = cj, the discontinuity distorts locally the ridge. 
We first consider m = for simplicity. Using the following development of the 
error function, 

w{z) = 1 + ^z -z^- —^z^ + 17 + • • • , 
we obtain the following expansion for the spectrogram: 



P5/(0, ^) = 2(1 + cosa;r) + x^ sina;T - 



TT 



1 + 1 + - cosa;T 
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sina;T' + 4a; 
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37r 
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where we have dropped the overall constant (7-^/2. The equation that deter- 
mines the extrema of the spectrogram at = is the following one: 



sin ujT — X 



2/2 
1 - - + 1 + - ) cosa;r 

TT V TT 



x~—= sina;^ 



+ 2x' 



l- — +\l + — ]cosujT 



+ ••• = 



(B.l) 



From a solution xq of this equation, we obtain an approximation of the ridge: 



V2 

U Xq. 

a 



One of the main features of a solution x^y of this equation is that it does not 
depend on how large is the time interval. It only depends on the fractional part 
of a;T/(27r). Moreover this error decreases like 1/cr as a is increased. Again we 
notice that when loT — mod 27r the ridge is located at ^ = a;. 
We neglect the quadratic and cubic term of Eq. (B.l) which is valid if ^ = loT 
mod 27r is close to zero. There is only one ridge at = which is approximately 
located at 



1 



Li) — 



sin LoT. 



a\/27r 

The ridge is above ^ — lo when Q is negative, and below ^ — lo when Q is 
positive. 

When = n mod 2n, there is an extremum at Xq = 0. However this solution 
is a minimum of the spectrogram. By considering the quartic term in the 
equation, we see that there are two other local minimum given by: 



2a' 



Both ridges are local maxima of the spectrogram. 

When the quadratic term of the expansion of the spectrogram vanishes (for 
coT Ri 1.7947 mod 27r), Eq. (B.l) shows that there are two extrema for the 
spectrogram: 



± 




Only ^_ is a local maximum of the spectrogram and is located below ^ — oo. 
We compute next an approximation of the ridge when u is increased. Using 
the expansion oi w, we obtain 



Pshu, 
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The equation that determines the ridge is the following one: 




yTT L TT V TT/ J L VT \ TT/ 



5 r 8 / 8 \ " 

-x'^{l-y^)^smu;T+2x^(l+y^) 1 - — + 1 + — cosu;T +--- = 




The approximation of the ridge is given by 



1 




sin LoT. 



The ridge decreases (or increases) with u if the ridge is above (resp. below) 
the main frequency uj. The time interval where the effect of the discontinuity 
at = is visible is approximately [0, r] where 



The length of this interval is independent of the total time length T and is 
proportional to a. If we increase a, the difference between the ridge and the 
frequency uj decreases as l/c but the time interval where this effect is visible 
increases as a. Figure B.l displays the ridge plots in three cases: The first panel 
is for ^ = where a flat ridge is expected, the second panel is for ^ = 7r/4 and 
(7 = 1, and the third panel is for 9 — tt/A and a — 1.5. 

For ujT = TT mod 27r, the extremum a; = is a minimum for y small and 
becomes a maximum for y > 1. Thus the typical time of the effect of the 
discontinuity is again r = a\/2. For y < I, there are two ridges approximately 
located at : 



These two ridges are symmetric with respect to ^(u) = u. Figure B.2 dis- 
plays the two ridges in that case. The asymmetry on this figure is due to the 
correction of higher order terms in x in the expansion of the spectrogram. 

B.2 Wavelet transform 

The effect of a discontinuity at m = due to a periodisation is approximately 
the same as for the windowed Fourier transform. The computations are very 
similar. The wavelet transform of / is given by 



(j\/2. 




Wf{u,s) = {a^n) 




Again two terms are important in the sum n — —1 and n — 0. 
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Jo 

Next we assume that T is large and we express the wavelet transform using 
the complex error function w. 



2 \ 1/4 



Wf{u,s)'^ \ ^ ^ w{x + iy) + e^'^ w{—x — iy) 

where x — a{suj — ri)/y/2 and y — —u/ {say/2). The main difference with 
the windowed Fourier case is that now y depends on the frequency (through 
the scale). For at] sufficiently large and s sufficiently close to rj/cu, we get an 
approximation of x and y: 



where ^ = r]/s. Thus the results of the windowed Fourier transform can be 
transposed for the wavelet transform by replacing a by ar]/uj. 
Using the expansion of the error function, we obtain an approximation of the 
value of the maximum of the wavelet coefficients and hence an approximation 
of the ridge. For u — Q and for uT ^ mod 27r, this value is : 

a; 

^ — uj -1= sm ujT. 

(Tr/v27r 

In the case of the wavelet ridge, we notice that the difference is proportional to 
the frequency, i.e. the the error in determining the frequency u will be smaller 
for low frequencies than for large ones. 
For u close to zero, we get an approximation of the ridge: 

f[u) — (jj 1= 1 — ^ „ „ sm Lol . 

The typical width in time of the error in determining the frequency uj is given 

by 

£777^2 

r a; . 

LO 

We notice that increasing the parameters of the wavelet decrease the amplitude 
of the error but increases the typical length time where this error is visible. 
The wavelet method cannot determine frequencies such that r > T/2, i.e. 
frequencies such that 

a;<2V^^ 
T 
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Fig. B.l. Wavelet ridge of a function /(t) = e*"^* with uj = 1/2 on an interval [0,T] 
for (a) uT = 207r, (5) = 207r + 7r/4 and a = 1, and (c) = 207r + 7r/4 and 
(7 = 1.5. 

In summary, the fact that the trajectory is known on a finite time interval 
creates in general a discontinuity on the boundaries of this interval. This dis- 
torts the ridge on a small time interval. It creates one or two ridges; when 
there are two ridges, it means that the main original ridge is a minimum of 
the spectrogram. When there is one ridge, it can be below or above the main 
ridge depending on the value of the fractional part of ujT/{2t^). We have seen 
that the amplitude of this distortion is inversely proportional to the size of 
the window of the Fourier transform and that the typical time interval during 
which this distortion is visible is proportional to this size of the window. 
Figure B.3 represents the influence of the discontinuity in the time-frequency 
plane. The gray region is the portion of the time-frequency plane that has to 
be discarded due to the influence of the discontinuity on the boundaries. For 
the windowed Fourier transform, this error does not depend on the frequency 
whereas for the wavelet transfrom, the distorsion influences greatly the low 
frequencies. 
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Fig. B.2. Wavelet ridge of a function f{t) = e*'^* with to = 1/2 on an interval [0, T] 
for LoT = 207r + tt. 
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Fig. B.3. Schematic representation of the finite time effect on the computation 
of the ridges: The gray region is the part of the time-frequency plane where the 
discontinuity of the signal between t = and t = T influences the computation of 
the instantaneous frequencies (a) for the wavelet transform and (6) for the windowed 
Fourier transform. 
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